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Abstract 

150 y63.rS after the discovery of quaternions, Hamilton's conjecture that 
quaternions are a fundamental language for physics is reevaluated and shown 
to be essentially correct, provided one admits complex numbers in both 
classical and quantum physics, and accepts carrying along the intricacies 
of the relativistic formalism. Examples are given in classical dynamics, 
electrodynamics, and quantum theory. Lanczos's, Einstein's, and Petiau's 
generalizations of Dirac's equation are shown to be very naturally formulated 
with biquaternions. The discussion of spin, isospin, and mass quantization is 
greatly facilitated. Compared with other formalisms, biquaternions have the 
advantage of giving compact but at the same time explicit formulas which 
are directly usable for algebraic or numerical calculations. 
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Notations 



In general scalars, vectors and quaternions are represented by the following types : 

• Scalars or biscalars (elements of M or C) : 
Preferably lower case Roman or Greek: a, 6, c, a, p... 

• Vectors or bivectors (elements of M.^ or ) : 



• Quaternions or biquatemions (elements of H or B) : 
Preferably upper case roman or Greek: A, C, A, E... 

Square brackets are used to emphasize that bracketed quantities (e.g., [x] for a 
variable, or [ ] for an operand within an expression) are quaternions, or to represent 
quaternions as [scalar; vector] pairs : 

• [x]^[s;v\^ s + v e H or B 

Angle brackets are used to emphasize that bracketed quantities, e.g., {S), are 
scalars, or to restrict quaternions to their scalar part : 

9 {s + v) = s eMorC 

Please note the use of the following operators : 

• or ()~ quaternion conjugation (bar or minus) 

• 0* imaginary conjugation (star) 

• O"*" biconjugation (plus) 

• ( )~ order reversal or ordinal conjugation (tilde) : (AB)^ — B^A^ 



Ao B scalar part of quaternion product AB : Ao B = (AB) = §>[AB] 
AaB vector part of quaternion product AB : AaB^ Y[AB] 



Any arrowed character: a, B, V, u, n, ... 



( )~ transposition (transpose) : 




( ) proper time derivative (dot) : () = 7-p — () 
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With c the velocity of light (we use generally the convention c = 1) the most 
common four-dimensional physical quaternions are : 

• four-position : X — [ict;x] 

• four-gradient : V = [^; ^] = [di^t, V] 

• four- velocity : t/ = 7[1; —i(3] 

• energy-momentum : P — [E; —icp\ 

• generalized momentum : 11 = [H; —icn] 

• electromagnetic four-potential : A = [ip; —iA] 

• charge-current density : J — [p; —ij] 

• electromagnetic induction bivector : F —[0;E + iB], F~ = [0; -B — iB] 

• electromagnetic excitation bivector : G = [0;D + iH], = [0; D — iH] 

• electromagnetic induction tensor : F{) = {)F + ) 

• electromagnetic excitation tensor : G{) — {)G -\- 

Lorentz transformation quaternions operators are 'mathcal' characters : 

• Rotation-spinor : TZ{) =exp(|^a)[] = [cos(|^); sin(|^)a][ ] 

• Boost-spinor : B() = cxp(i|?/6)[] = [cosh(|7/);isinh(|7/)6][ ] 

• Lorentz-spinor : C{) — BR[ ] 

Quatemionic and quantum mechanical scalar products : 

• Angle "bra" and "ket" delimiters are used to emphasize the range of symbols 
over which a scalar part is calculated : 

{A\B) = {AB) = ^[AB] 

• Formal quantum mechanical Hilbert spaces 's scalar products are emphasized 
by using double "bra-ket" symbols, i.e., ((...)) instead of (...). E.g.: 

where ( )"•" is biconjugation, not Hermitian conjugation denoted by ( )^. 
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1 Introduction 



For a contemporary physicist, the name Hamilton is primarily associated with 
what is known as the "Hamiltonian" formulation of dynamics. Although math- 
ematically equivalent to other formulations such as the Lagrangian formalism, 
Hamilton's method provides a description of a classical system which has the con- 
siderable advantage that the problem can easily be "quantized," i.e., generalized 
from classical to quantum physics. For this reason, in all text books, the operator 
H appearing on the right-hand side of the Schrodinger wave equation 

z^l-^ = H^, (1) 
ot 

is called the "Hamilton operator" or, simply, the "Hamiltonian." 

On its own, this formulation of dynamics discovered by Hamilton in 1834 
(which allows problems of optics and problems of mechanics to be worked out 
interchangeably) is enough to put Hamilton among the greatest physicists of all 
times, at least at the level of Newton or Maxwell, and not very far from Einstein. 
At this Conference, however, we celebrate another major discovery of Hamilton: 
The invention of quaternions in 1843. 

But quaternions were not Hamilton's only important discovery in algebra: 
Complex numbers were first. Indeed, in 1835, Hamilton had already found a 
mathematically appealing and consistent way of interpreting the so-called "imagi- 
nary numbers." By considering pairs of ordinary numbers, and defining a suitable 
multiplication rule, he showed that all operations that could be made with ordinary 
numbers could also be made with his number doublets. 

However, as it often happens in the mysterious process that we call "discovery," 
Hamilton had a peculiar mental image in his mind when he was thinking of algebra: 
While he imagined geometry of a science of "space," he conceived algebra as the 
science of "pure time" [1 J, and therefore understanding imaginary numbers meant 
for him coming closer to understanding the essence of time. Thus, when Hamilton 
was thinking of algebra, his mental image was that of a physicist, an image of 
somebody whose ambition is to discover the laws of inanimate nature and motion. 

Therefore, when after many unsuccessful attempts Hamilton finally succeeded 
in generalizing complex numbers to quaternions (which require for their repre- 
sentation not just two but four ordinary numbers) he definitely believed to have 
made a very important discovery. This conviction, however, Hamilton would not 
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include in his scientific writings. But in his correspondence^ and in his poetry]^ 
Hamilton made it plain that he really thought he had discovered some synthe- 
sis of three-dimensional space, the vector-part, and time, the scalar-part of the 
quaternion. 

The astonishing fact is that indeed quaternions do foreshadow "our four- 
dimensional world, in which space and time are united into a single entity, the 
space-time world of Einstein's Relativity," Ref. [2], page 136. In effect, as sci- 
ence advances, more and more evidence accumulates, showing that essentially all 
fundamental physics results can easily and comprehensibly be expressed in the 
language of quaternions. If that is so, then the often-made criticism that Hamilton 
had "exaggerated views on the importance of quaternions," Ref. [2], page 140, was 
ill-founded. And therefore, contrary to what is often said, Hamilton was right to 
have have spent the last twenty-two years of his life studying all possible aspects 
of quaternions. 

However, what Hamilton did not know, and could not have known at his time, 
is that quaternions would only become really useful in applied and theoretical 
physics when problems are dealt with in which relativistic and quantum effects 
play an essential role. In such problems, in effect, it is not so much the "real 
quaternions" that Hamilton was mostly studying which are useful, but the so-called 
biquaternions, which are obtained when every four components of the quadruplet 
are allowed to become complex numbers. This is not to say that Hamilton's work 
on real quaternions was vain. Quite the contrary: Most algebraic properties of real 
quaternions that Hamilton so carefully studied can be carried over to biquaternions . 

But, for the more practical purposes of non-relativistic or non-quantum physics 
and engineering, it is true that operations with quaternions are not sufficiently 
flexible. We are therefore in a fortunate position today, that after the simplification 
of the somewhat cumbersome notations used by Hamilton, we have at our disposal 
the modern vector notation introduced by J.W. Gibbs. Using this notation, we 
can now work with quaternions much more easily than with Hamilton's original 
notation!^ In particular, either we separate the quaternion into its scalar and vector 
parts, separating or mixing freely "scalars" and "vectors" [3], or use it as a whole, 
especially when we deal with the fundamental aspects of theoretical physics, as in 
this paper. 

'"My letter relates to a certain synthesis of the notions of Time and Space, ..." Cited in Ref. IT], 
page 149. 

^"And how the One of Time, of Space the Three, Might in the Chain of Symbol, girdled be." 
Cited in Ref. [l ^page 192. 

''When using the modern vector notation with quaternions it is important to keep in mind 
Gibbs's redefinition of the sign of the scalar product: a • & = — (ab) — ~aob = — §[a&] . 



6 



Hence, if we now anticipate what we will develop in the sequel, and jump 
to our conclusion, we are indeed going to show that while quaternions are not 
the panacea for solving all possible physical or mathematical problems, they do 
nevertheless provide an extraordinary ly powerful framework for any problem in 
which some four-dimensional or quantal aspect of our physical world intervenes. 

That this is so, and why it is so, is mysterious. As Wigner stressed in his 
often quoted article entitled "The Unreasonable Effectiveness of Mathematics in 
the Natural Sciences" [4J, the biggest mystery is, possibly, the fact that once a 
particularly efficient mathematical scheme has been found for the description of 
some often crude physical experiment, it turns out that the same mathematical 
tool can be used to give an amazingly accurate description of a large class of 
phenomena. This is what we have discovered: Once the biquaternion formalism 
is taken seriously as a language for expressing fundamental physical laws, it so 
happens that more and more phenomena can be predicted by simply generalizing 
the accepted results while staying within the frameworks of quaternions. 

To do so, one has however to follow a few simple guide-lines. For instance, 
the time variable should always be written as a pure imaginary number 
and, consequently, derivation with respect to time should always be written d/dit. 
Hence, the fundamental space-time variable X and the corresponding four-gradient 
V will always be written as 

X=[ict;x], (2) 

A second rule is that there should be no hidden "i;" in other words, that the 
imaginary unit "z" should always be explicit, and that imaginary conjugation 
should always apply to all "f's. This means that contrary to the convention 
of some physicists, e.g., EISl, one should not use the so-called "Hermitian-" 
or "Pauli-units," but only the real quaternion units defined by Hamilton, which 
together with the scalar "1" have the advantage to form a closed four-element 
group, which is not the case with the "Pauli-units." 

If these two rules are followed, one discovers that there is one and only one 
"f ' in physics and in mathematics; that imaginary conjugation can always be 
given a consistent interpretation (either in classical or quantum physics); and that 
while "i'H is necessary in quantum theory, "i" is also a very useful symbol in 
classical physics because it often contributes to distinguish quantities which are of 
a different physical nature. 



'^Or "ici," where c is the velocity of light, to make explicit that the time and space variables 
have different physical dimensions. 

^Or, more precisely, a "complex structure" ||6l. 
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To conclude this introduction, let us summarize our main point: If quaternions 
are used consistently in theoretical physics, we get a comprehensive and consistent 
description of the physical world, with relativistic and quantum effects easily taken 
into account. In other words, we claim that Hamilton's conjecture, the very idea 
which motivated more then half of his professional life, i.e., the concept that 
somehow quaternions are a fundamental building block of the physical universe, 
appears to be essentially correct in the light of contemporary knowledge. 



2 Some definitions and properties of 
quaternion automorphisms 

Let us review some definitions and properties of the four basic quaternion linear 
automorphisms, i.e., the non-trivial involutions of H, the field of quaternions, or 
B the algebra of biquatemions. 

The first two are quaternion conjugation, which reverses the sign of the vector 
part, and imaginary conjugation, which replaces the scalar and vector parts by 
their imaginary conjugate 

Q^Q=[s--% (4) 

Q^Q* = [s*-v*]. (5) 

Quite often in practice ( ) and ( )* are used in combination. Following Hamilton's 
usage of the prefix bi- we call this third involution biconjugatioi^ and use for it 
the symbol ( )+ 

Q ^ Q+ = [s*; -V*] = (Q)*. (6) 

In the same spirit we call a complex vector a bivector (rather than a "six- vector"), 
but we will refrain from using the term "biscalar" suggested by Hamilton for 
complex numbers. 

Using these three involutions we have the following definitions: 

• Q is a scalar if Q = Q • Q is a vector if Q = —Q 

• Q is real if Q* = Q • Q is imaginary if Q* = —Q 

• Q is bireal if = Q • Q is antibireal if = —Q 
^Rather than "Hermitian conjugation," symbol ( Y, as it is often improperly called. 
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When operating on a quaternion expression, quaternion conjugation reverses 
the order of the factors. Thus 

AB = BA and = (7) 

The last non-trivial involution, order reversal (or ordinal conjugation), is more 
subtle and requires some explanations for which it is best to return to Hamilton's 
"Elements of Quaternions," and more specifically to a note added by C.J. Joly in 
1898, at the end of section ten [7J, Vol.1, p. 162. 

Starting from the set of quadruplets of real or complex numbers, the quaternion 
algebra is obtained by requiring their product to be associative, and the division 
to be feasible always, except possibly in some singular cases. Then, writing 
two quadruplets A and B as scalar- vector doublets [a; a] and [b;b], and using 
contemporary vector notations, their product has the following explicit form 

[a; a][b;^ = [ ab + p a ■ b ; ab + ab + q a X b]. (8) 

The two constants p and q are related by the equation 

q^ + p^ = 0, (9) 

which shows that there is some residual arbitrariness when defining the product 
of two quadruplets. For instance, taking p = —1, q can be equal to either +1 or 
— 1. On the other hand, taking p = +1, q may be +i or —i. Thus, the choice 
p = +1 corresponds to the Pauli algebra. But, as we have already said, we will keep 
Hamilton's choice, p = —1, which is also more fundamental because it corresponds 
to the Euclidian metric in the case of real quaternions, and to Minkowski's metric 
in the case where the four-dimensional space-time position vector is written as 
in formula ©. Moreover, with Hamilton's choice, the imaginary conjugate of a 
product is equal to the product of the imaginary conjugate of the factors: This 
dispenses of the special rules which are sometimes necessary when using "Pauli 
units.'izi 

In short, the arbitrariness in the sign of p is connected with the signature of 
the metric, and the choice of the metric determines the sign of p and the usage of 
"z" when defining physically meaningful four-dimensional quantities. In effect, 
since the square of the norm of a quaternion A is by definition its product by its 
conjugated quaternion, we have |v4p = A A = + p\v\'^ . Therefore, for a given 
signature, the choice of the sign of p is immaterial because one can always multiply 
the vector part of all quaternions by "i" in order to get the desired signature. 

^Such special rules are also necessary in the standard "7-matrices" formulation of Dirac's 
theory because they are based on the Pauli matrices which contain an algebraic ^(—1) which 
should not be mistaken with the "i" of the complex scalars that multiply the "7"s. 
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The arbitrariness in the sign of q is due to the non-commutativity of the 
quaternion product. Indeed, changing the order of the factors A and B is equivalent 
to changing the sign of q. The involution associated with the changing of this sign 
is called "order reversal" (or simply "reversal") and is designated by the symbol 
( When biquaternions are used to represent physical quantities in space-time, 
since q is the sign associated with the vector product, there is a close connection 
between order reversal and space inversion. However, contrary to the case of 
p, there appears to be no invariant overall criterion to decide for the sign of g@ 
Therefore, in accordance with the principle of relativity, one has to make sure that 
fundamental physical entities are "order-reversal covariant" (or simply "ordinal 
covariant"), i.e., that they do not arbitrarily depend on the sign of q. 

In this respect a last point is of importance: Whereas the problem of signature 
is common to all formalisms, order reversal is specific to quaternions and Clifford 
numbers^ and therefore something that should be carefully considered when any 
Clifford algebra is used in physics. For this reason, because reversal was not 
properly considered — for instance, in defining fundamental quantities as "ordinal 
invariant" — many authors using quaternions in physics have met with problems^ 
Indeed, as we will see, reversal plays an essential role when writing fundamental 
equations of physics. 



3 Irreducible representations of biquaternions: 
Spinors and bispinors. 

Before going to physical applications, let us remind (without proof) a few important 
elementary theorems concerning the irreducible decompositions of biquaternions: 

Theorem 1 Any normed real quaternion, i.e., R = R* such that |-Rp = 1, can be 
written 

a ad 
R = exp(-a) = [cos(-); sin(-)a], (10) 

where a is a real number called the angle and a a unit vector called the axis. 

^The operator O"', which is postfixed to an operand, should not be confused with symbol 
denoting the "dual" of an entity, i.e., ( ). 

''This observation relates to the experimental fact that while time flows in only one direction, 
space is not oriented. 

'''"Order reversal" is related (but in general not identical) to "reversion," one of the three basic 
involutions defined on any Clifford algebra. 

"(Note added in 2002.) See, for example, Refs. ||62ll63l . 
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Theorem 2 Any normed bireal quaternion, i.e., B = B^ such that \B\'^ = 1, can 
be written 

B = exp(i-6) = [cosh(-);isinh(-)6], (11) 

where y is a real number called the rapidity and b a unit vector called the boost 
direction^ 

Theorem 3 (Main biquaternion decomposition theorem) Any biquaternion Q 
with non-zero norm, i.e., \Q\'^ ^ 0, can be written 

Q = cRB, (12) 

where c is a complex number, R a normed real quaternion, and B a normed bireal 
quaternion. 

Theorem 4 Any biquaternion S with zero norm, i.e., 15*^ = 0, in which case S is 
called singulaO can be written as a product of three factors 

Q = rRa, (13) 

where r is a real number, R a normed real quaternion, and a a primitive nullquat 

a = i[l;zz7], (14) 

where V is a real unit vector, which has the property of being an idempotent, i.e., 
a^ = a. 

Theorem 5 Multiplying a nullquat from one side by any biquaternion does not 
change its primitive nullquat. 

Since two real parameters are needed to fix the direction of a unit vector, we 
see from (flOl) and (fTT)) that three parameters are necessary to represent a normed 
real or bireal quaternion. Similarly, a general biquaternion (fT2)) requires eight 
parameters, while six suffice for a nullquat (fT3l) . 

Theorem[5]provides the basis for defining spinors in the quaternion formalism. 
In effect, for a given primitive nullquat a, the left- (or right-) ideal forms a four 
parameter group that is isomorph to the spin | spinor group [|9l|. As a consequence, 
we have the spin ^ decomposition theorem which establishes the link between 
Dirac's bispinors and biquaternions: 

^^A bireal biquaternion is called a minquat by Synge fS^ . 
'^A singular biquaternion is called a nullquathy Synge |8|. 
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Theorem 6 (Spin 1/2 decomposition theorem) Relative to a given primary nul- 
Iquat a, any normed biquaternion Q can be written as a bispinor, i.e., as the sum 
of two conjugated spinors 

Q = Ri(7 + R2a, (15) 
where Ri, R2, are two real quaternions. 



4 Relativity and quaternionic tensor calculus 

Usually, when special relativity is introduced, one does a lot of algebra in order to 
work out the somewhat complicated formulas of Lorentz transformations. Even 
when quaternions are used for this purpose, the proof of the equivalence of the 
quaternion formulas with the usual ones is rather complicated [lOJ. In the present 
paper, our intention is to develop the fundamental concepts and present the main 
results without giving the details of the proofs. 

We start therefore from the fundamental ideas of relativity, covariance, and 
tensor calculus which are that all observers are equivalent for writing the physical 
laws, and that all meaningful physical quantities should have well defined trans- 
formation properties when going from one observer to another one0 Hence, if 
Hamilton's conjecture is correct, i.e., if indeed biquatemions can be used as ele- 
mentary building blocs of theoretical physics, any meaningful physical quantity 
should be writable as a simple explicit quaternion expression which should have 
the same form in all referentialsl!^ In other words, the components which in or- 
dinary tensor calculus are represented by symbols such as tfj ;, where the various 
indices show how the physical quantity varies in a change of referential, should be 
replaced by quaternionic monomials QRST... where the indices are replaced by 
some convention making the variancei^ of each factor in the monomial explicit. 

Let us take, for example, the general Poincare transformation law. This is a 
change of referential which corresponds to the affine function 

Q' = AQB + C, (16) 

where (a priori) A, B and C are any kind of quaternion expressions. If A and B 
are functions of the four-position vector X, and the translation term C is zero, we 
have a local Lorentz transformation, and if A and B are independent of X a global 
Lorentz transformation. 

''^For an enlightening introduction to tensor calculus and its relation to quaternions see ifTTI . 
See also Ref. [2J page 140. 

^^I.e., "reference frames" (gallicism). 
^^I.e., "transformation law" (gallicism). 
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The most basic tensor quantity is obviously the four-vector such as, for exam- 
ple, the space-time vector X given in formula Q. Four- vectors have the following 
variance 

V = CVC^, (17) 
where C is restricted by Einstein-Minkowski's condition 

W = l/'y' = invariant scalar. (18) 

Using (fTTl) we see that this condition implies 

CC = C*C^ = 1. (19) 

Thus, the quaternion C which represents the most general Lorentz transformation, 
is simply a biquatemion of unit norm. By theorem IH such a transformation can 
therefore be decomposed into a product C = TZB where 7^ is a rotation and B 
a boost. To find the explicit form of the boost, we apply (flTl) to the velocity 
four-vector U. Then, transforming from the rest-frame (in which U = 1) to a 
moving frame, we find 

U' = BB+ = BB = ^[1; -i(3] , (20) 

where P = v/cis the relative velocity of the moving frame and 7 = (1 — /?^)~^/^ 
the Lorentz factor. We see therefore that the Lorentz boost is a kind of quatemionic 
square-root of the four-velocity. 

But, the four-vector is not the most simple non-trivial covariant quantity. 
Hence, for a spinor, there are four possible transformation laws (or eight, if one 
takes order reversal into account) 

= jOSi , 5*2 = , iSg = C 5*3 , = S^C^ . (21) 

In fact, as it is immediately seen, these four types of variances are the counterparts 
to the four basic variances of tensor/spinor calculus: Contra- or co- variance, 
dotted or undotted indices. But, here, the last three can be deduced from the first 
by means of the automorphisms (), ()*, and ()+. This leads directly to the general 
idea of quatemionic tensor calculus: Any time some covariant four-dimensional 
quantity is introduced, the only possible new variance it may obtain is the result of 
operating with one of the three basic involutions, possibly combined with reversal 
( )~. Multiplying these quantities, and alternating the variances by making use of 
quaternion conjugation, one obtains more complicated tensors. For example, 

f = V[l^l^2], (22) 
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is a six-vecto^ which has the variance T' = C*TC^ . In fact, since C*C^ = 
CC = 1, the scalar part {V1V2) of V1V2 is an invariant, while its vector part 
T = Vi A V2 is a complex vector such as, for example, the electromagnetic field 
bivector F = E + iB = VAA. 

Of course, when changing referential, besides the proper transformations such 
as (fT6l) . there are also the so-called improper ones which involve space- or time- 
reversal, complex conjugation, and order reversal. These can be taken care of in 
the quaternion formalism, using in particular the four basic involutions, so that, in 
the spirit of relativity, covariance can also be insured with respect to them. For 
example, in the case of tensors constructed by multiplying four-vectors of odd 
parities, a tensor that is ordinal invariant will also be of odd parity. 

In summary, as long as one remains within four-dimensional space-time and 
works with biquatemions, it is possible to achieve the same power and flexibility 
as with ordinary tensor/spinor calculus without having to manipulate explicitly a 
large number of indices. This is possible in the realms of both special and general 
relativity [[13. 

5 Classical dynamics and Hamilton's principle 

Classical mechanics is a domain in which Hamilton himself found many brilliant 
applications of real quaternions. Just think, for example, of his very elegant 
and general resolution of the Kepler problem, in which the quaternion formalism 
leads directly to the conservation of angular momentum and of the misnamed 
"Runge-Lenz" vector0 

Here, however, we will consider the classical dynamics of a system of point 
particles, without specifying a priori whether it is a relativistic or non-relativistic 
problem. The fundamental concept is then the "Hamiltonian," which is a scalar 
function of time, position x, and canonical conjugate momentum tF. This Hamil- 
tonian function H can be merged together with the canonical momentum into one 
bireal quaternion that we call the four-Hamiltonian 

11= [H{t,x,Tf);-iTf{t,x)]. (23) 

For example, for a system point particles in a time independent external field, 11 

'^Similarly to the concept of four-vector, the concept of six-vector refers to a variance, not just 
to the fact that such objects are necessarily bivectors which have six real components. 

''^See Ref. [7|, Vol II, art. 419, pages 298-299. The angular momentum and Runge-Lenz 
vectors are denoted by /3 and e, respectively. 
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is the sum 

n = ^P„ + g„A(XO, (24) 

where _P„ = [En, —icpn] is the energy-momentum of each particle, g„ their electric 
charge, and A = —iA], the external electromagnetic four-potential. 

The equations of dynamics can then be expressed in a number of equivalent 
forms, the simplest one being possibly the action postulate which states that the 
four-Hamiltonian derives from an invariant scalar function, the action S{t, x) : 

VS = m. (25) 

Since S" is a scalar, operating on both sides with V and taking the vector part we 
find _ 

VAn = 0. (26) 

This is, written in quaternions, the condition for dU to be a total differential, or, 
equivalently, for 11 to be an exact one-form, i.e., 

dXoU = 0, (27) 

which, for any given two fixed point Xi and X2, is the same as Hamilton 's principle 

5 I Z^on = 0. (28) 

Moreover, for a system in which H does not depend explicitly on time, (|27l) is also 
fully equivalent to Hamilton 's equations of motion: 

S=^H , ^ = -^H. (29) 

dTT ax 

The sequence of transformations we have gone through may look like a succession 
of trivialities. This, in fact, is not the case. Had we not put the "i" in front of 
7? in (|23l) . we would not have been able to work out these results. Moreover, 
while all expressions are formally covariant, they are the same whether we assume 
the kinematical expression for P to be relativistic or not. Hence, although we 
have done nothing more than rewriting well known results, we see that complex 
quaternions provide a compact and convenient framework for writing the equations 
of Hamiltonian dynamics, and that the resulting expressions are automatically in 
relativistic covariant form. It is the same with Maxwell's equations: As will be 
recalled in the next section, writing them down in compact quaternionic form 
requires the use biquaternions. 
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To show what happens if we introduce relativity suppose for example that we 
apply (l25l) to a single particle. The fourth component, i.e., time or energy, is no 
longer an independent variable. For instance, we have now the relativistic identity 
|Pp = and (l25l) can be rewritten as 

(iV^ + eA){iVS + eA) = ml (30) 

This is Hamilton- Jacobfs equation for a particle in an electromagnetic field. 



6 Maxwell's equations and quaternionic analyticity 

After the casting of Lorentz transformations into quaternion form, one of the first 
modern applications of biquatemions was the rewriting of Maxwell's equations in 
191 1 by Conway Ca, and in 1912 by Silberstein [[Ml, as 



VA = F, 
VF = -47rJ, 



(31) 



where F = E + iB is the electromagnetic field bivector, J = [p, —ij] the source 
charge-current density, and A = —iA] the electromagnetic potential for which 
the Lorenz gauge, i.e., V o A = 0, is assumed. This very compact form allows 
many calculations to be done very effectively. In particular, as shown by Silberstein 
in 1913, the energy-momentum tensor of the electromagnetic field, i.e.. Maxwell's 
stress-energy tensor, has a very simple explicit quaternionic form [fTSl 



47rT() = iF+[]F = -iF~[]F, (32) 

where the free space [ ] corresponds to the position of a quaternionic argument^ 
When this tensor is used to calculate the flow of energy and momentum through an 
hypersurface the result is automatically covariant, and there is no "4/3 problem" 
as with the obnoxious Poynting vector [16J. Moreover, using ([3T]) to calculate its 
divergence, one immediately obtains Lorentz's force-density equation 

dp 

= -T(V) = -\{JF + F-J) = -F{J), (33) 



^'This convention, due to Hamilton and promoted by Silberstein, Conway and Synge fSl, gen- 
eralizes Dirac's "bra-ket" notation to biquaternions. Its value stems from the speed of calculation 
which derives from the simplicity of the composition rule: a[ ]a' b[ ]b' = ab[ ]b'a' . Moreover, it 
provides a clear distinction between "numbers" (or "vectors") Q, and "functions" (or "operators") 

Q[]- 
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which shows that the quaternion form of the electromagnetic field tensor is 
-^0 — |([]-^ + -^~[])' ^ physical object which should not be confused with 
the electromagnetic field bivector F = E + iB, or its reverse F~ = E — iB, a 
non-trivial distinction first made in 1955 by Kilmister [fTTl . 

A most interesting idea suggested by Maxwell's equations in quaternion form 
was developed by Lanczos ifTSl in his PhD thesis of 1919. In effect, Maxwell's 
second equation in vacuum, VF = 0, is the direct generalization of the Cauchy- 
Riemann analyticity condition from two to four dimensions. It is therefore natural 
to envisage classical electrodynamics as a biquatemionic field theory in which point 
singularities are interpreted as electrons [fT9l . In this case the field at some point 
X is calculated by means of the appropriate generalization of Cauchy's formula 
in which the integration contour becomes an hypersurface surrounding the 
point 

where R = Y — X and |-Rp = RR. This generalization of complex analysis has 
been extensively studied by Fueter in the case of real quaternions [20] and more 
recently extended to biquaternions and higher dimensional Clifford algebras EH- 
This formalism can now very efficiently be applied to standard problems, such as 
the calculation of retarded potential and fields [22] . 

7 Spinors in kinematics and 
classical electrodynamics 

Spinors are increasingly often used in classical physics and relativity [|23ll . How- 
ever, possibly the first significant use of spinors in classical physics was made 
in 1941 by Paul Weiss [|24||. the particularly brilliant PhD student of Max Born. 
By "significant" we mean that Weiss's applications of quaternions was not just 
rewriting an otherwise known result in quaternion form. In fact, Weiss gave 
an independent interpretation and derivation of an important physical law: The 
Lorentz-Dirac equation. 

Weiss's starting point was the fact that the quaternion formalism provides 
explicit formulas which are difficult to obtain by the ordinary methods of analysis. 
For instance, in kinematics, taking the square root of the four-velocity as in (|20l) 
is the same as making the spinor decomposition of the four-velocity. An explicit 
formula for the four-acceleration is then obtained by taking the total proper-time 
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derivative on both sides of 

Z = U = BB+, (35) 

which leads to 

Z = U = -BAB+, (36) 

where, as shown by Weiss, the invariant real vector A is the acceleration in the 
instantaneous rest-frame. 

Similarly, since null-four-vectors can explicitly be formulated with biquater- 
nions, one has explicit formulas for the light-cone and retarded coordinates, i.e, 

X~Z = 2i^BaB+, (37) 

where ^ = —iU o {X — Z) is the invariant retarded distance from the position of 
the charge Z to the space-time point X, and a an idempotent such as ([14)) with u 
pointing from Z io X. 

In his paper, Weiss does not speak of spinors. On the contrary, he makes it 
clear that his decomposition has nothing to do with Dirac's bispinors. But what 
he does is exactly the kind of spinor decomposition we use today, e.g., in general 
relativity. 

Weiss's application is to show that in this formalism the flow of energy and 
momentum through a hypersurface surrounding a point charge in arbitrary motion 
can be calculated exactly using Silberstein's form of the Maxwell tensor (|32|) . 
He then proceeds to find the world-lines for which the energy-momentum flow is 
stationary, and discovers that the resulting equation of motion is nothing but the 
Lorentz-Dirac equation [i25l 

mc^U = lie^iuijU + U)- \e{UF + F^U). (38) 



8 Lanczos's generalization of Dirac's equation: 
Spin and isospin. 

About one year after Dirac discovered his relativistic wave-equation for spin ^ 
particles, Lanczos [26J published a series of three articles in which he showed 
how to derive Dirac's equation from the more fundamental coupled biquatemion 
systeni^ 

VA = mB, 
VB = mA. 



(39) 



2" We write m for mc/h taking c = fi, = 1 for simplicity. Note that Lanczos could have taken 
the reverse of (|39]| as his fundamental equation: AV — mi3~ , B^'V ~ mA. ((39l ^) 
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Obviously, Lanczos was inspired by his previous work with quaternions [fTSl. 
Indeed, comparing with (|3T]) . it is clear that (|39l) can be seen as Maxwell's equations 
with feedback, and that, following Lanczos ll27l . this feedback can be interpreted 
as a distinctive feature of massive particles. However, there is a problem. In 
Dirac's equation, the wave function is a four-component bispinor, while A and 
B are biquaternions with four complex components each. This is the "doubling" 
problem that puzzled Lanczos a lot, as well as others who later tried to cast Dirac's 
equation in quaternion form [|28l . 

The first step towards a Dirac bispinor is to postulate that A and B have 
spinor variances, i.e., that A' = CA and B' = C*B, which leaves the possibility 
of making a gauge transformation, i.e., a right-multiplication by some arbitrary 
biquatemion G. Then, to get Dirac's spin | field, i.e., a bispinor, Lanczos had to 
make the superposition (see theorem[6]in Sec. [3]) 

D = Aa + B*a. (40) 

Here a is an idempotent such as (fT4l) . where V is any unit quaternion. Compar- 
ing with (fTSl) . we see that o has the effect of projecting out half of A, which 
added to another half of the complex conjugated of B, gives a Lorentz covariant 
superposition that obeys the wave equation 

VD = mD*iv. (41) 

This equation, to be called the Dirac-Lanczos equation, is precisely equivalent to 
Dirac's equation. It will be rediscovered by many people, in particular by Giirsey 
[|29l and Hestenes [130]| . While equivalent to other possible forms, (|4TI) has the 
considerable didactic advantage of making "spin" explicit. Indeed, the unit vector 
on the right-hand side shows that Dirac's equation singles out an arbitrary but 
unique direction in ordinary space: The spin quantization axisEH 

Using this equation, it is easy to construct and study the various covariant 
quantities which are important in quantum electrodynamics. For example, the 
conserved probability current is J = DD^, and Tetrode's energy-moment tensor 
is 

T{) = {{[]V)D)tuD+ - Dti){D+{V[])) - {AM)DD+, (42) 

where is the potential of an external electromagnetic field. 

However, the superposition (|40l) is not the only one leading to a spin | field 
obeying equation (|4T]) . As shown by Giirsey in 1957, if (l40l) represents a proton, 

^'(Note added in 2006.) Even more important, the complex conjugation on the right-hand side 
exphcitly shows that the Dirac field is fermionic, rather than bosonic as the Maxwell and Proca 
fields rt271l63l. 
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the neutron is then OTl [32ll 



N = {Aa- B*(r)iT, 



(43) 



where r is any unit vector orthogonal to u. 

Hence, Lanczos's doubling is nothing but isospin. Giirsey's articles had a 
tremendous impact p7] and inspired ideas like chiral symmetry and the sigma 
model [132] . Indeed, "internal" symmetries such as isospin are explicit and trivial 
in Lanczos's double equation (l39l) . while only space-time symmetries are explicit 
in Dirac's traditional 4x4 matrix formulation or the biquatemionic formulation 
(|4TI) . Unfortunately, except in his PhD dissertation, Giirsey made no reference to 
Lanczos's work, and Lanczos never learned that he had anticipated isospin in 1929 
already! 

Since the fundamental fields are A and B, while D and are the physically 
observed ones, it is of interest to find the most general gauge transformations on A 
and B which are compatible with the superpositions (|40l) and (|43]) . In fact, these 
transformations form a group that was discovered in another context by Nishijima 
[|33l and which has the following explicit representation 



By direct calculation, one finds that while A and B transform under Gn, D 
transforms as exp(— z/a) and as exp(— z7/3), respectively, so that the system 
(|39l) describes two particles of equal mass but different electric charges, such 
as the proton and the neutron. Hence, by just trying to write Dirac's equation 
in quaternions, one is automatically led to discover the existence of isospin, a 
fundamental feature that indeed is found in nature^ 

9 Proca's equations and the absence of 
magnetic monopoles 

When we wrote Maxwell's equation (|3T]) we made the implicit assumptions A = 
A"", i.e., that A was a ordinal invariant fundamental four- vector. If we try now 
to put a mass term on the right of the second Maxwell equation, and therefore 
introduce a "feedback" to get the wave equation for a massive spin 1 field, we find 
that Maxwell's equations have necessarily to be generalized to the following form 



Gn = crexp(ia) + aexp{i(3). 



(44) 




(45) 



'For more details, see Sec. 7 of Ref. ll35l . 
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This is, written in biquaternions, the correct spin 1 wave equation discovered in 
1936 by Proca Il34l . As with Dirac's equation it is easy to write in quaternions the 
conserved current and the energy-momentum tensor 

2J = + + (...)~, (46) 

SttTO = F+[]F + mM+n^ + (...)~, (47) 

where {■■.)'" means that the expression has to be completed by adding the reverse 
of the part on the left. Hence, the current and the energy-momentum are bireal 
and ordinal invariant four- vectors, as it should be. 

Now, just as we derived Dirac's equation from Lanczos's equation (|39l ) by 
making the superposition (l40l) . Proca's equation (|45|) can also be derived from 
(|39l) by adding the second Lanczos equation to its reverse equation. If this is so, 
what then is the meaning of the equation obtained by subtracting Lanczos's second 
equation from its reverse, assuming that the potential A is ordinal invariant 

VF-F~V = ? (48) 

Obviously, this is just the part of Maxwell's equation which specifies that there are 
no magnetic monopoles! Hence, if Lanczos's system (|39l ) is taken as the funda- 
mental equation from which Dirac's and Proca's equations are derived. Maxwell's 
equation is obtained by taking the = limit in (|45l) . and (|48]) insures the 
absence of magnetic monopoles. 

However, if we would have assumed that A = —A^ instead of A = +A^, 
we would have found another fully covariant field equation, only differing from 
Proca's by the fact that a minus-sign would replace the plus-sign in (1451) : In fact, 
the correct equation for a massive pseudo-vector particle. Therefore, as shown by 
Giirsey in his PhD thesis [|29l , the wave-equations of all scalar and vector particles, 
and of all pseudo-scalar and pseudo-vector particles, are just degenerated cases of 
Lanczos's fundamental equation (|39l) . 

10 Einstein-Mayer: Electron-neutrino doublets 
in 1933! 

When he wrote his 1929 papers on Dirac's equation, Lanczos was with Einstein 
in Berlin. In 1933, Einstein and Mayer (using semi-vectors, a formalism allied 
to quaternions) derived a spin ^ field equation (in fact, a generalized form of 
Lanczos's equation) predicting that particles would come in doublets of different 
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masses [|35l[36]| . The idea was that the most general Lagrangian for quaternionic 
fields, to be called the Einstein-Mayer-Lanczos (EML) Lagrangian, should have 
the form 



(49) 



L = § A^VA + B+VB - {A^BE+ + B+AE) + (...)" 

where i? G B. The field equations are then 

VA = mBE+, \ 
VB = niAE. J 

which reduce to (|39l) when E = m. In the general case, the second order 
equations for A ox B become eigenvalue equations for the mass (the m factor 
appearing in the argument of exp(^im(Et — p ■ x)) of plane wave solutions. This 
generalization is obtained by the substitutions Am — > AE and Bm Bm in the 
Lagrangian leading to (l39l) . Therefore, mass-generation is linked to a maximally 
parity violating field. 

There are two basic conserved currents: The probability current J, and the 
barycharge current K 



J = AA+ + BB+ , K = AEA^ + BEB+. (51) 

Keeping E constant, J is invariant in any non-Abelian unitary SU{2) ® f/(l) 
gauge transformation of A or bE On the other hand, K is only invariant for 
Abelian gauge transformations which also commute with E, i.e., elements of the 
general Nishijima group (l44l) such that E and a commute. 

Of special interest are the cases in which E is also a global gauge field. The 
first such gauge is when E is idempotent. One solution of (l50l) is then massive 
and the other one massless: An electron-neutrino doublet! The most general 
local gauge transformations compatible with (ISTI) are then elements of the unitary 
Nishijima group f/Ar(l, C) combined with one non-Abelian gauge transformation 
which operates on either E and A, or E and B, exclusively. This leads directly to 
the Standard model of electro-weak interactions llJTll . 

The second fundamental case is when E is real: E = E*. The eigenvalue 
equation is degenerate and the two masses are equal. As shown by Giirsey 1321 . 
equation (l50l) describes a nucleonic field, non-locally coupled to a pseudoscalar 
sigma-pion field. 

Again, as in the examples shown in the previous sections, Hamilton's con- 
jecture seems to be realized. Einstein-Mayer's generalization (l50l) lifts the mass 
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In biquaternions, these transformations have the form G{ ) — [ Je*"^ exp{^9a) 
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degeneracy of Lanczos's original equation (|39l ) and leads to electron-neutrino dou- 
blets and weak interaction on one hand, and to proton-neutron doublets and strong 
interaction in the form of the well known charge-independent pion-nucleon theory, 
on the other hand [|27l@ 

11 Petiau waves and the mass spectrum of 
elementary particles 

One of the central problems of contemporary physics is the question of the origin 
of the mass of the elementary particles. As we have seen, by replacing the mass 
term in Lanczos's equation by some biquaternionic parameter, Einstein has been 
able to show that elementary particles come in doublets of different masses. In 
fact, in 1930 already, Lanczos wondered whether a theory in which the mass term 
is replaced by a variable would not simultaneously solve the problem of mass 
quantization and that of infinities in field theory [|39l . 

Well, as nobody after Einstein and Mayer seemed to have taken Lanczos's 
suggestion seriously, one had to wait until 1965 (about the time when Gell-Mann 
and Zweig proposed the idea of quarks) for somebody to reinvent the concept. This 
year, in complete independence from mainstream research, the French physicist 
Gerard Petiau wrote a system of equations which may precisely give a solution to 
the problem of the mass of the electrons and quarks [i40il . 

Although Petiau was thinking in very general terms, considering complicated 
couplings between particles of various intrinsic spin, it is easy with quaternions to 
write his fundamental equation in the case of spin | particles [f35l. It just amounts, 
in the spirit of Lanczos's feedback idea, to adding a third equation to Einstein's 
system (l50l) in order to close it: 

VA = BC, ] 

VB = AC, \ (52) 
VC = AB. J 

Here A and B are the usual Lanczos spin | fields, and the scalar C an additional 
Einstein-Mayer field of spin 0. Because the system is now closed, it becomes 

^'*(Note added in 2006.) Soon after the publication of Einstein and Mayer's papers, Pauli asked 
his assistant Valentine Bargmann to study them. Bargmann found that the solutions of Eq. dSOl l 
consist of doublets of particles with different electric charges and masses, including particles of 
zero mass ||38l . It is unfortunate that Pauli, who had postulated the existence of the neutrino in 
1927, did not understand the importance of Einstein and Mayer's papers, which contain all the key 
ingredients of the contemporary Standard Model if E is interpreted as the Higgs field 1621 . 
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non-linear, and the solutions are much more constrained than with any usual linear 
type of wave equations. 

For instance, the single-periodic de Broglie waves that quantum mechanics 
associates with a particle become double-periodic Petiau waves [41]. Instead 
of being linear combinations of sin (2;) and cos(z) functions, these waves are 
superpositions of elliptic functions sn(z,k), cn(z,k),etc. A very appealing feature 
of Petiau waves is that their dependence on the modulus interpolates between pure 
de Broglie waves (for A; = 0) and pure soli tonic waves (for A; = 1): A beautiful 
realization of the wave/particle duality of quantum mechanics. Moreover, both 
the amplitudes and the proper mass (the /i factor appearing in the argument of 
sn(^lj,{Et — p- x),k), for example) will be quantized. 

The most interesting thing, however, is happening when, in order to quantize 
the system, the Hamiltonian function is constructed. Taking, for example, A 
as the fundamental field Petiau showed that, in terms of the first integrals, the 
Hamiltonian has the very simple form [40] 

H = Co/A;2, (53) 

where k is the modulus of the elliptic function, ji the proper mass, and Co some 
constant. The exciting thing is that the Hamiltonian, and thus the total energy in 
the field (i.e., for a single particle, the effective mass) scales with the fourth power 
of /i. 

In effect, in 1979, Barut discovered a very good empirical formula for the 
mass of the leptons [42|. Assuming that a quantized self-energy of magnitude 
^a~^MeC^N'^, where = 0, 1, 2, is a new quantum number, be added to the 
rest-mass of a electron to get the next heavy lepton in the chain e, /x, r, Barut 
got the following expression (where a = 1/137) 

n=N 

M{N) = M,{l + la-^J2^^)- (54) 

n=0 

The agreement with the data of this rather simple formula is surprisingly good, the 
discrepancy being of order 10~^ for /i and 10^'^ for r, respectively [|43ll . In order 
to get the masses of the quarks 11441 . it is enough to take for the mass of the lightest 
quark M„ = Me/ 7. 47 . Again, as can be seen in Table 1 the agreement between 
the theoretical quark masses and the "observed" masses is quite good, especially 
for the three heavy quarks. 

Since we have just seen (|53l) that the energy of a Petiau field is scaling with the 
fourth power, we are inclined to think that there might indeed be a fundamental link 
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between such non-linear fields and the theory of the mass of quark and electrons. 
If this is so, what about the factor 7.47 ? 

There are two non-trivial exceptional cases for elliptic functions: The harmonic 
case, k = sin(|), and the equianharmonic case, k = sin(Y^). It is very plausible 
to associate the former with leptons, and the latter with quarks. Indeed, in either 
case, the corresponding elliptic functions exhibit several unique symmetry and 
scaling properties, which come from the fact that in the complex plane their poles 
form a modular aggregate with f or | symmetries. Since according to (52) the 
masse is proportional to k'^, the electron to quark mass ratio is then equal to 
[sin(f)/sin(fL)]2^7.47. 

But this is now very close to pure speculation, and in any case on the frontier 
of contemporary research [|44ll . Nevertheless, it is interesting to see how far, just 
following Hamilton's conjecture, one can go in the direction of a unified picture 
of fundamental physics. 

Table 1 
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Comparison of electron and quark masses in MeV/c^ calculated with Barut's 
formula (|54l) to measured masses from Ref. h43\l . (Note added in 1996: The 
observation of a sixth quark of mass in the range of 160 '000 to 190 '000 MeV/c^ 
has been reported at the beginning of 1995.) 

12 Quaternions and quantum mechanics 

Some of those who have been following us on this upwards trail, starting from 
Hamilton's principle and ending with a possible solution to the problem of quark 
masses, might be surprised that just a short section is dedicated to quaternions and 
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quantum mechanics. To these we say — in the spirit of Hamikon's particle- wave 
duality — that everything we have done can easily be recast in the jargon of "wave 
mechanics," so that, in this perspective, we have been doing quantum mechanics 
all along. 

In effect, the quantum predicate is no so much in the field equations we have 
been discussing in this paper, than in the interpretation of the field's amplitudes 
(the "wave functions:" complex number in the non-relativistic case, biquaternions 
in the relativistic one). As Feynman clearly stated in a review of the principles 
of quantum mechanics: "It has been found that all processes so far observed 
can be understood in terms of the following prescription: To every process there 
corresponds an amplitude; with proper normalization the probability of the process 
is equal to the absolute square of this amplitude" Ref. [|45 J . page 1. 

Take, for example, the Dirac-Lanczos equation ((4T)) . rewritten here in the case 
where there is an external electromagnetic field 

VD = {mD* + e'A^,^D)iU. (55) 

Comparing with (1), we see that the Hamiltonian is the following operator 

H{) = -V[]-m[]*iv-eA,^[]iV. (56) 

Then, following Feynman's prescription, we have to normalize the amplitude D. 
Since the probability current is the conserved four-vector DD^, a suitable norm is 

{{D\1\D)) = JJJd'V {D+[1]D) = 1, (57) 

where the dummy operator "1" may be replaced by the operator corresponding 
to the physical quantity whose expectation value is to be calculated. Hence once 
the above prescription has been accepted as a postulate there is not much mystery 
left in quantum theory, and it is straightforward, at least in principle, to give the 
quantum interpretation of the field equations presented in the previous sections]^ 

At this point, it is worth mentioning that in the past decade a kind of a silent rev- 
olution has been occurring with respect to quantum theory. Increasingly, quantum 
and classical theories are seen as part of the same theory: Barut shows that wave 
mechanics can be formulated without h [|46l . Lamb suggests that Newton could 
have invented wave mechanics [|47l . and several major investigations show that 
the whole apparatus of the so-called "second quantization" of fields is redundant 
fi48il . Indeed, as is trivially shown in the way Lanczos's equation ([391 ) generalizes 

^^It is remarkable that it is the use of complex conjugation in expressions such as the Hermitian 
product (ISTI i that distinguishes quantum theory from classical physics. 
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Maxwell's equation [1311, Maxwell's theory can be interpreted as a quantum theory 
without h. For instance, calculating by means of Silberstein's energy-momentum 
tensor (|32l) the electromagnetic field's energy-momentum density, or the Lorentz's 
force density (l33l) . is the same as applying the quantum rule (|57l) . Moreover, 
the Hamiltonian is simply the operator H( ) = |(V[ ] — [ ] V) just like in Good's 
quantum interpretation of Maxwell's theory ||49ll . 

What, is then the main contribution of biquatemions to quantum theory? Pos- 
sibly, the clear disentanglement of "i" and "/i," the two elements which have 
been traditionally associated with quantum mechanics. Indeed, looking at the 
Schrodinger equation ([T]), these two elements appear together as a combined fac- 
tor. In Lanczos equation [|39l . however, "i" appears in the scalar part of the 
four-gradient V, and h in combination with the mass on the right-hand side had 
we not used the convention h = 1. Thus, if Hamilton's conjecture is true, "z" is 
definitely associated with "time" (i.e., Hamilton's intuitive conception of imag- 
inary numbers) and is associated with "mass" or, more precisely, with the 
particle aspect of waves, i.e., lumps of energy localized in space 

It remains, in conclusion, to stress that the power of Hamilton's conjecture 
seems not to have shown its limits yet. By this we allude to the numerous 
investigations in "quaternionic quantum mechanics" which have occurred since 
the birth of wave mechanics. 

Indeed, not to mention the work of Lanczos 1*261 [39l. quaternionic and other 
more general algebraic generalization of quantum mechanics have been actively 
studied since 1928 already [|52l . The best known sequel to this work is possibly 
what Jordan ll53l started in 1932 and which culminated in the famous article [|54ll 
in which the first "Jordan algebra" was described. To give another example of the 
breadth of research in these early days, we mention that the theory of operators in 
quaternionic Hilbert spaces was the subject of a PhD thesis in 1935, and that the 
name "Wachs space" was proposed for such spaces [[55||. 

A new impetus was given in the 1960's, mainly after the work of the group 
around Finkelstein and Jauch [56] at CERN, followed by others [57], up to the 
synthesis soon to be published by Adler ll58l . All these developments contemplate 
the possibility that the complex numbers of contemporary quantum theory may 
have to be replaced by quaternions or biquatemions in some more fundamental 
theory. However, it may well be that Nature is satisfied with complex number as 
the fundamental scalar field, and that in this respect a single and commutative "i" 

^^(Note added in 2006.) The origin of "i" in quantum theory is explained by Lanczos as follows: 
"For reasons connected with the imaginary value of the fourth Minkowskian coordinate ict, the 
wave mechanical functions assume complex values" lISTI p. 268]. 
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is enough and playing some essential role that the current experimental situation 
seems to favor [|59l[60l . 



13 Conclusion 

The physics is in the mathematical structure, not in the, formalism: What are then 
the advantages of using a formalism such as Hamilton's biquaternions? 

• Biquaternions are whole symbols, i.e., they compound between one and eight 
real numbers which belong to a single (or a few related) tensor quantity(ies) 
so that many formulas written in biquaternions are simpler than their standard 
matrix, tensor, or higher rank Clifford numbers counterparts. In general, 
they enable to dispense of at least one level of tensor indices, and quite often 
to reduce a few indices tensor to a single entity. 

• Biquaternions are expressive, i.e., being elements of the simplest non-trivial 
Clifford algebra they provide neat and explicit formulas in final form, which 
are therefore directly amenable to symbolic or numerical calculation, with 
pencil and paper, or with a computer. 

• Biquaternion formulas are suggestive, i.e., they often indicate the correct 
way of generalizing a result, or how to relate seemingly independent results. 

• Biquaternions provide a unifying formalism, e.g., they enable a fully con- 
sistent use of complex numbers in both classical and quantum physics; they 
lead to expressions that are very similar in both Galilean and Lorentzian 
relativity; they are very effective in formulating the physics of the current 
"Standard model" of fundamental interactions [|27l [35l[62l ; etc. 

However, some words of caution are in order: Since biquaternions are whole 
symbols it is important to take care of the problems specifically associated with 
such symbols. For instance, in order that biquaternion expression have a well 
defined tensor character they have to be constructed from elementary biquater- 
nions that have such a character. Moreover, special care is required because of 
noncommutativity and of the need for biquaternion expression to be ordinal co- 
variant. For example, the truly correct form of Maxwell's equations is not the 
Conway-Silberstein expression (|3TI) . but the gauge and ordinal invariant system 



V A A = F, 
|(VF + F-V) = -47rJ, 



(58) 
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which like Proca's equation (|45l) does not require the supplementary conditions 
(Vv4) = 0, A = y4~, 3 = 3"", i.e., the Lorenz gauge and the requirements that A 
and J are ordinal invariant biquatemions. It is is precisely because such problems 
were not properly understood and cured that biquatemions failed, at the beginning 
of the twentieth century, to become a widespread language for physics. 
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